A new possible quark-hadron mixed phase in protoneutron stars 
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The phase transition from hadronic matter to quark matter at high density might be a strong 
first order phase transition in presence of a large surface tension between the two phases. While 
this implies a constant-pressure mixed phase for cold and catalyzed matter this is not the case 
for the hot and lepton rich matter formed in a protoneutron star. We show that it is possible to 
obtain a mixed phase with non-constant pressure by considering the global conservation of lepton 
number during the stage of neutrino trapping. In turn, it allows for the appearance of a new kind of 
mixed phase as long as neutrinos are trapped and its gradual disappearance during deleptonization. 
This new mixed phase, being composed by two electric neutral phases, does not develop a Coulomb 
lattice and it is formed only by spherical structures, drops and bubbles, which can have macroscopic 
sizes. The disappearance of the mixed phase at the end of deleptonization might lead to a delayed 
collapse of the star into a more compact configuration containing a core of pure quark phase. In 
this scenario, a significant emission of neutrinos and, possibly, gravitational waves are expected. 
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The possibility of the formation of a quark hadron 
mixed phase in neutron stars, especially during the first 
stages after their birth, has been widely discussed in the 
literature and interesting associated signatures were pro- 
posed [n,a i in. In Ref. [|| a delayed formation of the 

quark phase was found, after the deleptonization of the 
star, which can trigger the collapse of the protoneutron 
star to a black hole. An observation of a supernova neu- 
trino signal with an abrupt cessation of the signal would 
be a clear confirmation of this scenario. Recently, a sim- 
ilar result was obtained in Ref. Q where a supernova 
simulation for a 100 Mq progenitor star was performed. 
The appearance of quark matter reduces the delay be- 
tween the core bounce and the collapse to a black hole. 
In both above mentioned paper the MIT bag model was 
used to describe quark matter and a large value of the 
bag constant was adopted resulting in a large critical den- 
sity for the appearance of quark matter. By assuming 
a small bag constant value, which still allows to obtain 
maximum masses for hybrid stars compatible with obser- 
vations, it was found in Ref. [8| that quark matter can 
form already during the early post bounce phase of a 
core-collapse supernova. Interestingly, the phase transi- 
tion leads to the formation of a shock wave propagating 
outwards which provides, in some cases, the explosion 
of the supernova even in spherical symmetry. Moreover 
a new neutrino burst is emitted with a typical time de- 
lay with respect to the first neutronization burst of a 
few hundreds of milliseconds. In all these studies the 
Gibbs construction [9( is adopted to describe the quark 
hadron mixed phase but finite size effects are not included 
i.e. one assumes that the surface tension a between the 
hadronic phase and the quark phase is vanishingly small. 
Presently, the value of a is not known and the possibil- 
ity that it has a lar ge v alue, of say a ~ 100 MeV/fm 2 , 
cannot be excluded [lOfl . If this is indeed the case, the 
mixed phase equation of state must be computed in- 
cluding the effects from the surface energy eg and the 
Coulomb energy ec of the finite geometrical structures, 
which are dubbed the "pasta phases" [ll|, . The op- 



timal size and shape of a structure at a fixed density 
is determined by the minimization of the energy which 
gives the well known relation: £5 = 2ec- It was pointed 
out in Refs. [H,[Ii| that also the effect of charge screen- 
ing and the rearrangement of charged particles in pres- 
ence of the Coulomb interactions must be taken into ac- 
count for a realistic description of the mixed phase. In 
fact local charge neutrality is almost recovered at lengths 
much larger than the Debye screening length. There- 
fore the mixed phase window shrinks considerably and it 
app roaches the constant-pressure Maxwell construction 
EH • This effect in turn would imply the absence of 
the mixed phase in cold and deleptonized hybrid stars. 
A calculation of finite size and charge screening effects 
in the mixed phase for protoneutron star matter, which 
has fixed entropy per baryon and fixed lepton fraction, 
has not yet been performed. Here we want to consider 
this possibility but instead of solving the system of Pois- 
son equations, which is necessary to obtain the density 
profiles of charged particles, we assume local charge neu- 
trality for the two phases as in the Maxwell construc- 
tion (this simple treatment is justified if a is sufficiently 
large). Consequently, we introduce two distinct charge 
chemical potentials: ^ for the hadronic phase and n q c 
for the quark phase and one baryonic chemical poten- 
tial fj,B corresponding to the global conservation of the 
baryon number. In a newly born neutron star however, 
besides the conservation of baryon number also the lep- 
ton number is conserved as long as neutrinos are trapped, 
i.e. during the first ~ 10 seconds of the life of the star. 
Neutrino trapping allows to introduce a chemical poten- 
tial fiL, associated with the additional conserved leptonic 
charge, which coincides with the chemical potential of 
neutrinos fij, = fx v (see [l5| for a general discussion on 
mixed phases). Let us consider the interface between 
the two phases: a charge separated interface is formed 
with a size of the order of the Debye screening length, 
~ 10 fm, with a layer of positively charged, electron de- 
pleted, hadronic matter on one side and a layer of quark 
matter with an excess of the electron on the other side 
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FIG. 1: The equations of state are shown for the case of EoS2 
for protoneutron star matter (dashed line) and for cold and 
catalyzed star matter (solid line). The dots indicate the onset 
and the end of the mixed phase in both cases. 



FIG. 2: Density profiles for a star with a baryon mass of 
1.9M0 for the protoneutron star stage and the cold config- 
uration. The dots mark the onset and the end of the phase 
transition. 



(as discussed in [l(| for the CFL phase). The interface 
is stabilized by the resulting electric field. Notice that 
neutrinos, being not affected by the electric field, can 
freely stream across the interface. Consequently, lepton 
number is conserved only globally. This additional glob- 
ally conserved quantum number has similar effects as the 
global charge neutrality condition adopted to model the 
phase transition for vanishing values of the surface ten- 
sion. We consider here the "standard" conditions of a 
newly born neutron star [3j: the matter has a fixed lep- 
ton fraction Yj, = (n e + n v )jnB = 0.4 and fixed entropy 
per baryon S/Nb — 1 where n e , n u and tib are the elec- 
tron, neutrino and baryon number densities and S is the 
entropy. The conditions of chemical-equilibrium in the 
two phases read: 

Mn = M-B, M P = Ms + Mc (!) 
Mu = (Mb + 2/-40/3, fi d = (mb - Mc) A Ms = Md (2) 
Me = ML " MCj Me = ML - Mc> A = = l l L (3) 

where Mi(* = n ,P, e, v, u, d, s) are the chemical potentials 
of neutrons, protons, electrons, neutrinos and up down 
and strange quarks, respectively. The Gibbs conditions 
for mechanical and thermal equilibrium read: 

P h (HB,^, ML, T) = P q (HB,^ c , ML, T) (4) 

n h c - n h e = , n q c - n\ = (5) 

{l-X){n h & +n v )+x{nl + n v ) = Y L nB (6) 

(1 - X)s h + X-s q = S/Nn B (7) 

where P a , s a with a — h,q are the pressure and the 
entropy density in the hadronic and quark phases, x 
is the volume fraction of the quark phase and Uq are 
the electric charge densities of hadrons and quarks. As 
usual the baryon density is defined as follows: ub — 
(1 — x) n B + X n B- Notice that the eqs. (5) impose local 
charge neutrality of the two phases and eq. (6) global 
lepton number conservation. To calculate the equations 
of state of hadronic matter and quark matter we adopt 
the relativistic mean field model with the parametriza- 
tion TM1 for the former and the MIT bag model 



including perturbative corrections for the latter [ljj, ll8| . 
We set the masses of up and down quarks to zero and 
the mass of the strange quark to 100 MeV. We fix the 
constant which simulates the QCD perturbative correc- 
tions c = 0.2, and we select two values of the effective 
bag constant J5 e g in order to have a critical density for 
the phase transition in protoneutron star matter of ~ no 
and ~ 3no (where no = 0.16 fm~ 3 is the nuclear satu- 
ration density), corresponding to B^ 4 — 155 MeV and 

bH^ = 170 MeV. The two equations of state are labeled 
as EoSl and EoS2 for the two choices of the effective bag 
constant. In Fig. 1 we show the equations of state for 
matter in a protoneutron star (indicated with PNS) and 
for cold and catalyzed matter (indicated with "cold"). 
The remarkable result is that within the mixed phase 
constructed by solving eqs. (4)-(7), the pressure increases 
as a function of the density and a large range of density 
is occupied by the mixed phase. During deleptonization 
the pressure in the mixed phase gradually flattens and 
finally for deleptonized and cold matter one finds the 
usual result of a Maxwell construction with a constant 
pressure from the onset to the end of the phase transi- 
tion. We use now the above presented equations of state 
to study the structures of protoneutron stars and cold 
stars. In Fig. 2 we show the density profile for a pro- 
toneutron star and the corresponding cold configuration 
(assuming total baryon number conservation during the 
cooling and deleptonization of the newly born star.) The 
mixed phase, initially present in a ~ 2 km large layer of 
the star, gradually shrinks during the deleptonization of 
the star and finally disappears in the cold configuration. 
As a result, a sharp interface separating hadronic mat- 
ter from quark matter is obtained with a sizable jump of 
the baryon density. In Fig. 3 we show the mass-radius 
relations for the different cases. The black and grey (or- 
ange online) thick dashed lines correspond to hybrid pro- 
toneutron stars (EoSl and EoS2 respectively), the black 
and grey (orange online) thick solid lines correspond to 
the cold configurations. Neutron stars mass-radius rela- 
tions are also shown for comparison (thin curves labeled 
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FIG. 3: The mass-radius relations are shown for the different 
equations of state, EoSl, EoS2 and TM1 and for protoneutron 
stars and cold stars configurations (colors online) . The crosses 
correspond to the maximum mass configurations. 

with TM1). The new mixed phase appears in a pro- 
toneutron star because the pressure increases with the 
density. Therefore, we obtain stellar configurations with 
a core of pure quark matter, a layer of mixed phase and 
a layer/crust of hadronic matter for the case of EoSl and 
hybrid stars with only a core of mixed phase in the case 
of EoS2. The mixed phase cannot appear anymore in 
the star for cold and catalyzed matter because the pres- 
sure is constant and only configurations with pure phases 
are obtained. As discussed before, a sizable jump of the 
density occurs at the interface separating the two pure 
phases which affects the stability of the stars: at the 
onset of the phase transition the stars are gravitation- 
ally unstable and only if a sizable volume of the star is 
occupied by the quark phase the stars are stable. The 
mass-radius relations in this case correspond to the so 
called "third family" solutions, see 0, |20, [H HI for a 
detailed discussion of the properties of these stars. We re- 
mark that the new mixed phase introduced here has quite 
different properties with respect to the globally electric 
neutral Gibbs mixed phase: since the two phases are lo- 
cally charge neutral no Coulomb lattice with charged fi- 
nite structures of the two phases can form. Instead an 
amorphous phase with only spherical pasta structures is 
present (1-D and 2-D structures can indeed form only 
due to Coulomb interactions). Moreover the charge neu- 
tral structures have macroscopic sizes contrary to a glob- 
ally electric neutral mixed phase where the optimal size 
of the structures is limited by the Coulomb energy. To 
reduce the total surface energy the charge neutral struc- 
tures start to merge and asymptotically a full separation 
of the two phases will be obtained. Therefore in such 
a mixed phase no coherent scattering of neutrinos with 
pasta structures can take place [23| , as the neutrino wave- 
length is much smaller than the size of the structures. A 
detailed simulation of neutrino transport within this new 
mixed phase would be extremely interesting for the pos- 
sible implications on the neutrino signal of the changes of 
the structure of the star during deleptonization. Also the 
motion and the interactions of the drops/bubbles within 
the mixed phase, in presence of turbulence, might rep- 
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FIG. 4: The total baryon number of compact stars as a func- 
tion of the gravitational mass for EoS2 at different stages 
of the temporal evolution: protoneutron star (PNS, Yl=0A, 
S/N = 1), an intermediate stage with Y L = 0.25, S/N = 1 
and the cold configurations (cold). In the insert we show a 
magnification of the third family branch (solid line) and we 
also show the result for neutron stars (dotted line). HP, MP 
and QP denote the hadronic phase, the mixed phase and the 
quark phase respectively. 

resent an interesting source of gravitational waves [24| . 
Concerning phenomenology, an interesting possibility is 
a delayed transition of a protoneutron star in a third fam- 
ily star during/after deleptonization, which is outlined in 
Fig. 4. The plot shows the baryon number of the stars 
as a function of the gravitational mass for protoneutron 
stars (dashed line) and cold stars (solid line). For the 
sake of discussion an intermediate configuration is also 
included (dotted-dashed line) corresponding to partially 
deleptonized matter with Yj, = 0.25. The insert shows a 
magnification of the third family branch of cold hybrid 
stars. The letter A in the plot denotes the configuration 
of a cold neutron star which is unstable with respect to 
the collapse to a third family star, indicated by the letter 
B, with the same baryon number. The energy released in 
such a collapse (the difference between the gravitational 
masses of the two configurations at fixed baryon num- 
ber) is of the order of 10 ergs similar to values found in 
Ref. [25| . The letter C marks the maximum mass of cold 
hybrid stars. The letter D stands for the configuration 
of a protoneutron star for which the central density is 
equal to the density of the onset of the phase transition 
and a core of mixed phase is formed; the core of mixed 
phase increases with the central density and therefore 
with the mass of the stars. Stellar configurations with 
a baryon number lower than A are always composed of 
purely hadronic matter during the evolution of the star. 
Since C is smaller than D, for stellar configurations hav- 
ing baryon number between the labels A and C the corre- 
sponding protoneutron stars do not have quark matter in 
the core (neither pure phase nor mixed phase) but during 
the deleptonization, since the onset of the phase transi- 
tion decreases, at a certain point a core of mixed phase 
forms (e.g. at point E). As the deleptonization proceeds 
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the mixed phase gradually shrinks, a pure quark phase 
core starts to form and finally for fully dclcptonizcd mat- 
ter the mixed phase disappears and an hybrid star with 
pure phases is obtained. Depending on the detailed dy- 
namics of the formation of the pure quark phase core 
and the disappearing of the mixed phase it is possible 
that the evolution towards the final cold hybrid star con- 
figuration, for stars having a baryon number close to A, 
proceeds through a gravitational collapse (similarly to 
the transition from A to B). In that case the gravita- 
tional potential energy is released in a short amount of 
time and a burst of neutrinos and gamma rays can be pro- 
duced as proposed in Ref. [25j for the collapse to third 
family stars. In such a fast dynamics also gravitational 
waves might be emitted if non-radial modes are excited. 
On the other hand it is also possible that the evolution 
of the star proceeds through hydrostatical equilibrium 
configurations, most probably for stars having a baryon 
number close to C, and the gravitational potential energy 
is released gradually. No strong signature is expected in 
this case unless finite size effects do play an important 
role for the nucleation of the new phase and the hadronic 
star can be in a metastable state before converting into 
a hybrid star [26| . Finally, stars having baryon number 
larger than C will develop a core of mixed phase during 
deleptonization and will collapse to a black hole after the 
full deleptonization (similarly to the results of Ref. [6j]). 

In this paper we have shown that also in presence of a 
large value of the surface tension between quark matter 
and hadronic matter it is still possible to form a quark- 
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